Parachute Models Used in the Mars Science Laboratory 
Entry, Descent, and Landing Simulation 


Juan R. Cruz, 1 David W. Way, 2 Jeremy D. Shidner, 3 Jody L. Davis, 4 Richard W. Powell 5 
NASA Langley Research Center, Hampton, Virginia, 23681 

Devin M. Kipp, 6 Douglas S. Adams, 7 Anita Sengupta 8 
Jet Propulsion Laboratory, Pasadena, California, 91109 

and 

A1 Witkowski, 9 Mike Kandis 10 
Pioneer Aerospace, South Windsor, Connecticut, 06074 


An end-to-end simulation of the Mars Science Laboratory (MSL) entry, descent, and 
landing (EDL) sequence was created at the NASA Langley Research Center using the 
Program to Optimize Simulated Trajectories II (POST2). This simulation is capable of 
providing numerous MSL system and flight software responses, including Monte 
Carlo-derived statistics of these responses. The MSL POST2 simulation includes models of 
EDL system elements, including those related to the parachute system. Among these there 
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landing engines, and offloading. These models were kept as simple as possible, considering 
the overall objectives of the simulation. The main purpose of this paper is to describe these 
parachute system models to the extent necessary to understand how they work and some of 
their limitations. A list of lessons learned during the development of the models and 
simulation is provided. Future improvements to the parachute system models are proposed. 
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I. Introduction 

A N end-to-end simulation of the Mars Science Laboratory (MSL) entry, descent, and landing (EDL) sequence 1 
was created at the NASA Langley Research Center using the Program to Optimize Simulated Trajectories II 
(POST2). 2,3 This simulation is capable of providing numerous MSL system and flight software responses, including 
Monte Carlo-derived statistics of these responses. The MSL POST2 simulation (referred henceforth as 
“simulation”) provides the highest fidelity environment for assessing flight dynamic performance of the MSL 
spacecraft, and is the primary venue for verification and validation of flight dynamics requirements. The simulation 
includes numerous models of the various systems, including those related to the parachute system. This paper 
describes these parachute system models. 

The primary purposes of the parachute models within the simulation are to determine: 

• parachute opening forces, 

• aeroshell trajectory (e.g., velocity and altitude vs. time), 

• aeroshell dynamics (e.g., translational accelerations, rotational rates and accelerations), and 
• forces imposed by the parachute on the aeroshell. 

Parachute models were kept as simple as possible while still fulfilling their primary purposes as listed above. 
Thus, not all aspects of the parachute, its operation, and behavior are modeled in detail. For example, the trajectory 
and dynamics of the parachute deployment bag from mortar fire to suspension lines stretch are not modeled. Models 
were tuned to the available experimental data. Many of the models are empirical, and make use of data from 
previous similar missions to Mars (e.g., the Viking and Phoenix landers), and MSL-specific data. Thus, although 
the approach used in these models may be adapted to other missions, their use and the values of numerical constants 
in them should be examined critically before use. 

In a companion paper, 3 4 a reconstruction of the MSL parachute performance on Mars is presented and compared 
to some of the parachute models and results from the simulation. The descriptions of the MSL parachute system and 
its operation, as presented herein, are also included in reference 4. 
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II. Parachute Decelerator System 


A. System Description 

The Mars Science Laboratory used a single, supersonically-deployed disk-gap-band (DGB) parachute with a 
nominal diameter (Z> 0 ) of 21.35 m (see figures 1 and 2). This DGB parachute was similar in geometric proportions 

to those used by the Viking missions to Mars, 4 5 but 33 percent larger in nominal diameter, making it the largest 
parachute that has been deployed on Mars. 

The disk was fabricated from two types of fabric: a ripstop polyester of nominal 1.4 oz/yd 2 areal weight 

manufactured to a Pioneer Aerospace specification on the crown, and a ripstop nylon of nominal 1 . 1 oz/yd 2 areal 
weight per Parachute Industry Association specifications (PIA-C-7020B Type I for the white nylon fabric, 
PIA-C-7020C, Type I for the orange nylon fabric - see figure 2). All fabric on the disk was placed in the bias 
orientation. The various fabrics can be clearly seen in figure 2. The band was fabricated from the same ripstop 
nylon fabrics described above placed in the block orientation. The suspension lines were fabricated from 
Technora T221 braided cord. All other textile structural elements were fabricated from Kevlar 29. The confluence 
fitting (located at the single riser to triple bridle confluence point) was fabricated using titanium TI-6AL-4V plate 
and A286 steel bolts. 

Deployment was effected by a mortar 6 with a muzzle velocity of approximately 39 m/s. The parachute was 
qualified to be deployed at Mach numbers up to 2.3, with a not-to-exceed peak parachute aerodynamic opening 
force (design flight limit load, F dfll ) of 289 kN (65,000 lb). 7 

B. Concept of Operation 

The MSL EDL concept of operation, starting from just prior to parachute deployment, is shown in figure 3. 
Nominal operation proceeds as follows: 

• During entry the MSL aeroshell is balanced such that its center of mass is offset from its axis of symmetry. 
This causes the aeroshell to trim at a total angle of attack of approximately 16°. This nonzero total angle of 
attack trim angle produces lift, which allows MSL’s entry to be guided by banking the aeroshell. However, for 
parachute deployment, a near-zero total angle of attack is desired. Just prior to parachute deployment a series 
of six balance masses are ejected from the aeroshell at two-second intervals. The ejection of these balance 
masses brings the center of mass of the aeroshell to a location near its axis of symmetry, causing the aeroshell 
to fly at a total angle of attack of nominally less than five degrees - an acceptable value for supersonic 
parachute deployment. Deploying the parachute at a low angle of attack also reduces the aeroshell dynamic 
motions and minimizes side forces on the aeroshell during inflation. Concurrently with the ejection of the 
balance masses, the reaction control system (RCS) rolls the aeroshell to the bank angle required for later radar 
operation purposes. The ejection of the balance masses, and the rolling of the aeroshell, is known as the 
Straighten Up and Fly Right (SUFR) maneuver. During the SUFR maneuver the RCS is used to damp pitch 
and yaw oscillations. 

• After the SUFR maneuver is completed, the RCS is inhibited, the mortar is fired, and the parachute is 
deployed. Mortar fire is controlled by a velocity magnitude trigger. 8 

• Parachute deployment and inflation occurs during a time interval of approximately 2 s after mortar fire. 

• Ten seconds after mortar fire the RCS is re-enabled to damp aeroshell oscillations if necessary. 

• Fifteen to 30 seconds after mortar fire the heatshield separation conditions are met. Heatshield separation is 
controlled by a “dot product” velocity trigger 8 set to ensure an aeroshell Mach number less than 0.8 at 
separation. The RCS is inhibited just prior to heatshield separation. 

• Three seconds after heatshield separation the RCS is re-enabled to damp aeroshell oscillations if necessary. 

• The aeroshell descends under parachute for 40 to 180 s until the descent stage separation conditions are met. 
Descent stage separation is commanded by an attitude/velocity trigger. 8 The RCS is inhibited for a final time 
immediately preceding descent stage separation. 

• The descent stage free falls unpowered for 1 s, followed by another 0.2 s for main landing engines (MFE) 
warm-up, and an additional 0.8 s for de-tumbling and turning into the initial attitude for powered flight. One 
of the powered descent stage’s initial objectives is to perform a divert maneuver to avoid both short- and 
long-term re-contact with the backshell and parachute. The rover is placed on the surface of Mars by the 
powered descent stage operating in the sky crane mode. 

In figure 4 most of the parachute models described in this paper (left hand side) are shown vs. the EDL events 
(right hand side). The timeline is shown vertically from top to bottom. All times in figure 4 are approximate. 
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III. Modeling 


A. Basic Modeling Assumptions 

The simulation models the aeroshell and parachute as two rigid bodies as shown in figure 5. The parachute and 
single riser is modeled as one rigid body with its own set of mass properties (section III.D) and aerodynamic 
coefficients (section HIT). This rigid parachute/single riser body (referred henceforth as “parachute”) is connected 
to the aeroshell by a triple bridle. The single riser and triple bridle connect at the confluence fitting. This 
confluence fitting is not modeled explicitly as an independent mass element - its presence in this model is just 
geometric. The legs of the triple bridle are modeled as massless, tension-only, linear spring-dampers (section III.C). 
The parachute’s wind angles (i.e., angle of attack, angle of sideslip, and/or equivalents) and dynamic pressure are 
calculated at the parachute’s center of mass (see section V.B for recommendations regarding improvements in the 
calculation of these wind angles). The Mach number used for the parachute models, however, is that of the 
aeroshell’ s center of mass. This is done to eliminate unrealistic parachute Mach number oscillations arising from 
interactions between the elastic triple bridle and the area oscillations model. See section III.G for more details 
related to this Mach number replacement. Values for several of the system geometric parameters are presented in 
section III.B. 


B. Geometric Parameters 

Key geometric parameters of the parachute system are shown in table 1 and figure 6. Note that the nominal 
(reference) parachute diameter, D 0 , has an uncertainty associated with it. This uncertainty is used to randomly 
perturb the specific value of D 0 used for a given Monte Carlo trial in the simulation. Also note that the nominal 
parachute diameter, with perturbation, is used to calculate the nominal (reference) parachute area, S 0 . It is these 
perturbed values of D 0 and S 0 that are used in the simulation for calculations involving the trajectory, flight 
dynamics, and the parachute aerodynamic force at first full inflation, F fi . An additional, unperturbed set, of 
nominal parachute diameter and area, (^ 0 ) Des and (S 0 ) D ^, , are used for the calculation of the design parachute 
aerodynamic force at first full inflation, F fi Des , only, which is then used for comparison against the fixed design 
flight limit load, F dfll . The difference between F fi and F fi Des is discussed in section III.F. As-built geometric 
parameters of the MSL parachute are provided in table 2. 


C. Elastic Properties 

The only elastic elements in the parachute model are the triple bridle legs; the riser and parachute are assumed 
to be a single rigid unit. The legs of the triple bridle are modeled as massless, tension-only, linear spring-dampers. 
In this model the force on each triple bridle leg is calculated from: 


■k A T + r 

■ ik tbl lu^ tbl t ^ tbl 


d(^L rBL ) 


dt 


( 1 ) 


If F sd > 0 then F TBL = F sd else F TBL = 0 endif (2) 

Where, 

F sd is the force generated by a spring-damper system, 

k TBL is the spring constant (each leg), 

c TBL is the coefficient of viscous damping (each leg), 

A L tbl is the change in length of a triple bridle leg, 
t is time, and 

F tbl is the force on a triple bridle leg. 

The modeled elastic properties for each triple bridle leg are: 
k TBL = 1-712 x 10 7 N/m 

c tbl = 10.900 N # s/m 

The coefficient of viscous damping represents one percent structural damping at 5 Hz. The unstretched length of the 
triple bridle legs is given in table 1. These elastic properties become active at suspension lines stretch and are 
constant through the simulation. 
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D. Mass Properties 

As mentioned in section III. A and shown in figure 5, the parachute downstream of the confluence fitting is 
modeled as a single rigid body. The nominal mass properties of the parachute were calculated using a 
three-dimensional computer aided design (CAD) model using the as-measured masses of the components 
comprising the parachute, and an assumed inflated shape. This CAD model is axisymmetric. Thus, the parachute’s 
center of mass is on its axis of symmetry, the mass products of inertia perpendicular to the axis of symmetry are 
identical, and the mass products of inertia are zero. The nominal mass properties are given in table 3. For 
simulation Monte Carlo trials the mass properties were perturbed. The total perturbations were estimated using both 
correlated (i.e., perturbations due to assumed changes in the parachute inflated geometry) and uncorrelated 
components. A complete description of the parachute mass properties perturbations and how they were calculated is 
beyond the scope of this paper. The mass properties of the parachute become active at suspension lines stretch. 
There is no ramp-up of the mass properties during inflation. 

A portion of the parachute system mass properties (e.g., triple bridle, sabot) remains accounted for in the mass 
properties of the aeroshell. This portion of the parachute system mass properties is a small fraction of the aeroshell 
mass properties. Because of this, perturbations of their values are not necessary for Monte Carlo trials. 

E. Mortar Recoil Force 

The time at which the mortar recoil force starts is denoted as the mortar firing time for the purpose of this model. 
The mortar recoil force is imparted on the aeroshell as the parachute pack is ejected by the mortar at high speed. 
This mortar recoil force is modeled as a time-varying force of short duration. The mortar recoil force acts along the 
aeroshell axis of symmetry in the flight direction. No random perturbation of the mortar recoil force is included in 
this model. This model was created using data provided by General Dynamics Ordnance and Tactical Systems 
obtained from the MSL mortar qualification tests. 


F. Deployment, Inflation, and Opening Forces 

In this section the models for deployment, inflation, and opening forces are described. Because these models 
are closely interrelated, they are described in the same section. The models described in this section are empirical, 
and derived principally from Viking, 9 " 13 Phoenix, 14 and MSL test data. They are specific to MSL. Among the 
assumptions made in their development are that the parachute is deployed by a mortar while in supersonic flight, and 
that inflation is of the infinite-mass type. 

Deployment and Inflation 

For the purposes of this document, deployment and inflation are defined as follows: 

Deployment. The events occurring from mortar firing at time t = t MF , to suspension lines stretch at time t = t LS , 
that is t MF <t < t LS . For modeling purposes, the time at mortar firing, t MF , is defined as the time at which the 


mortar recoil force starts. Also for modeling purposes, the time at suspension lines stretch, t LS , is defined as 


the time at which the skirt of the canopy starts to emerge from the deployment bag. Except for the mortar 
recoil force, the parachute is not explicitly modeled during the time interval from mortar firing up to (but not 
including) suspension lines stretch. The time at suspension lines stretch, in seconds, is calculated from the 
equation 


H tb + Lr + L s 


( 3 ) 


The time at mortar fire, t MF , must be entered in seconds in equation (3). The length (H tb +L R +L s ^j is that of 
the stretched-out parachute at suspension lines stretch, and must be entered in meters in equation (3). The 
numerical value for this length can be determined from the data in tables 1 and 2. The parameter V MF _ LS is 
the average deployment bag speed from mortar fire to suspension lines stretch, and must be entered in meters 
per second in equation (3). Its value depends on the mortar muzzle velocity, V Muzzle , as given by 

Vhf-u, = W [0-875 + (0.060 -t/[0,l)- 0.030)] (4)* 


* The quantity l/[0,l) is a random variable with uniform distribution in the interval between 0 (inclusive) and 1 (not inclusive). 
Unless specifically called for, the value of l/[0,l) is not reused. 
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The muzzle velocity is given by 

V Muxte = 3 9 m/s + [4 m/s*t/[0,l)-2 m/s] (5) 

Note that V MF _ LS and V Muzzle are perturbed (random) parameters; they are calculated once for a given Monte 
Carlo trial. Equations (4) and (5) include both the nominal values (i.e., 0.875 in equation (4) and 39 m/s in 
equation (5)), and the perturbations (i.e., the quantities in parenthesis involving the t/[0,l) term) to be used 

in Monte Carlo trials. The calculation of the time from mortar fire to suspension lines stretch described by 
equations (3) to (5) is empirical and highly simplified. The accuracy of this calculation could be improved 
- see section V.A. 

At suspension lines stretch, the parachute “appears” in the simulation and its aerodynamic and inertial 
characteristics become active. This activation of the parachute occurs with the following initial conditions at 

^ = As • 

• The parachute’s airspeed vector at its center of mass, \ P oo , is 

V P ,.=V AS ,„ (6) 

where V A5;00 is the aeroshell’s airspeed vector at its center of mass. 

• The parachute appears at a perturbed (random) total angle of attack, a p Tot , and total angle of attack 
clock angle, (j) P Tot determined from 

a P , TM -T[ 0°, 8°, 12°] (7) f 

(pp,Tot = 360°* t/[0,l)-180° (8) 

Note that a p Pot and (j) P Tot are calculated once for a given simulation. Equations (7) and (8) include the 


perturbations to be used in Monte Carlo trials. In calculating \ Poo , a p Tot , and (p P Tot per equations (6), 

(7), and (8), it is assumed that the wind at the parachute’s center of mass is the same as that at the 
aeroshell’s center of mass. This is not exactly correct, but is accepted as a modeling approximation. If 
there is a wind gradient there can be a small difference in the wind at the two locations. In the case of 
MSL the high airspeed at parachute deployment, and the relatively small distance between the parachute 
and aeroshell centers of mass, makes this approximation reasonable. 

• The parachute body rotation rate vector is zero: Q chute = 0 . 

• The parachute is placed in a position consistent with its wind angles (i.e., a p Tot and (/) P Tot as defined 

above), the requirement that its axis of symmetry pass through the center of the triangle defined by the 
triple bridle attachment points on the aeroshell, and the requirement that none of the triple bridle legs be 
stretched beyond their unloaded length. Note that this positioning allows one or two of the triple bridle 
legs to be slack. 

The simulation implementation of the parachute activation model achieves the same objective as described 
above, but does so in terms of different variables. 

The triple bridle elastic properties (section III.C) and mass properties (III.D) become active (i.e., without 
ramp-up during inflation) at t = t LS . The aerodynamic properties (sections III.H and III.I) become active at 


t = t LS , and ramp-up as described in the opening forces model described below. 

Inflation. The events occurring from suspension lines stretch, t LS , up to, but not including, the first full 
inflation, t FI , of the canopy. That is, t LS < t < t FI . The time at first full inflation, in seconds, is calculated 
from the equation 


A 

V, 


inf 


( 9 ) 


where, 

V inf = 29.9 m/s + [7.0 m/s»t/[0,l)-3.5 m/s] 


( 10 ) 


t The quantity 7'[0°, 8°, 12°] is a random variable with triangular distribution whose lower limit is 0°, its mode is 8°, and its 
upper limit is 12°. Unless specifically called for, the value of T is not reused. 
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The time at suspension lines stretch, t LS , must be entered in seconds in equation (9). The nominal parachute 
diameter, D 0 , must be entered in meters in equation (9). Note that D 0 must include the perturbation defined 
in table 1. The parameter V inf is the average inflation speed, and must be entered in meters per second in 
equation (9). Equation (10) includes both the nominal value of V inf (i.e., 29.9 m/s), and the perturbation (i.e., 
the quantity in parenthesis involving the £/[0,l) term) used to calculate the value of V inf for Monte Carlo 
trials. 

Opening Forces During Inflation 

The forces and moments generated by the parachute during deployment are a small fraction of those 
encountered during inflation. Because of this, parachute deployment forces (except for the mortar recoil force) 
and moments are assumed to be zero in the present model. Parachute aerodynamic forces and moments during 
inflation are controlled by a time-dependent nondimensional variable, Kj , that is used as a multiplier for the 

parachute tangential, normal, and pitching moment coefficients (C r , C N , and C mo , respectively) . This 
multiplier is used during the time interval from t LS to t FI as shown in figure 7 and defined by 


For t LS < t < t FI (inflation): Kj = C x 


t — tj 




\ t FI * LS ) 


C x =l .407 (11) 

rij =4 

Equation (11) describes a prescribed area growth model where C x is the opening load factor and rij is the 
exponent of the inflation prescribed area growth power function. The value of C x is dependent on the drag 
coefficient model; thus, it is valid only for the drag coefficient and static aerodynamic coefficients model 
described in sections III.H and III. I. Note that for t > t FI the calculation of aerodynamic forces generated by the 

parachute proceed in the usual manner without using . 

Peak Opening Forces 

In this model the first full inflation is assumed to be coincident with the inflation first peak force. The 
parachute aerodynamic force at first full inflation, F PI is calculated from 


F FI =q a S 0 C x (C 2 T+ C 2 N f ( 12 ) 

where the values of the dynamic pressure, q x , C T , and C N are those at t = t FI . Note that the values of C T and 
C N in equation (12) need to incorporate the perturbations described in section III.I, be determined at the 
appropriate total angle of attack, and have incorporated the Mach Efficiency Factor described in section III.H. 
For comparison to the design flight limit load ( F dfll , a fixed value), he design parachute aerodynamic force at 

first full inflation, F fi Des , occurs at first full inflation, t = t FI , and is calculated from the equation 

F FI ,ues ={l + M De ^(S 0 ) Des C x {C D X igh ( 13 ) 

where (S 0 ) D ^ is an upper bound estimate of the nominal area (see table 1), and h is the upper bound of 

the drag coefficient model (section III.H). The values of q x and [C D ^ ^ f to be used in equation (13) are those 

at t = t FI . The value of {C Dq ^ f is determined using the Mach number at t = t FI . The leading term in equation 

(13), (1 + , is an additional multiplier for the calculation of F fi Des to add conservatism to the calculation 
of the design peak opening force. For the final set of simulations conducted in support of the MSF EDF, the 
value of A K Des was set to zero. The value of F fi Des is intended to be a conservative estimate of the parachute 

aerodynamic force at first full inflation on the parachute. Note that it uses conservative values of the nominal 
area (i.e., (S 0 ) Des as compared to S 0 ) and the drag coefficient (i.e., (C Dq ) h . h )- The MSF parachute was 


t This value of F n can be exceeded for t > t FI due to area oscillations. See section III.G. 
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qualified by test to a design flight limit load, F dfll , of 289 kN (65,000 lb). 7 For the results of a given set of 
simulation Monte Carlo trials (i.e., a set of entry conditions targeting a specific landing location) to be 
considered acceptable, the 99 th percentile value of F F1Des could not exceed F dfll . 

G. Area Oscillations 

Previous flight experience with DGB parachutes indicates that at Mach numbers ( M ) above 1 .4, large parachute 
force oscillations occur after the first inflation peak force (see, for example, reference 11). These force oscillations 
are at least partially related to changes in the projected area of the parachute, as the parachute undergoes collapse 
and re-inflation cycles, and are commonly referred to as area oscillations. A significant contributor to these area 
oscillations is interaction between the aeroshell wake and parachute flow fields. Because the MSL parachute was to 
deploy its parachute at a Mach number greater than 1 .4, and perhaps spend a significant amount of time in this flight 
regime (sufficient for several area oscillations), it was deemed necessary to include a parachute area oscillations 
model in the simulation. Including such a model was thought to improve the fidelity of the aeroshell dynamics 
results obtained from the simulation. 

The area oscillations model described here is synthetic, in the sense that it simulates area oscillations without 
modeling the underlying physics (which are poorly understood). It was, however, tuned to the available data from 
the Viking BLDT AV-4 flight test 11 and the Phoenix Mars flight. 14 " 16 The area oscillations model calculates a Mach 
number- and time-dependent nondimensional parameter, K AOM (M, t ) , that is used as a multiplier for the parachute 
tangential, normal, and pitching moment coefficients (C r , C N , and C mo , respectively). This multiplier is used 
during the time interval from first full inflation, t FI , up to the time at which the Mach number is 1.4. The value of 
K aom was calculated by a function consisting of a randomized sinusoidal oscillatory term constrained as follows: 
by an exponential function to drive K A0M to 1.0 as the Mach number approaches 1.4, by limiting the upper value of 
K aom , and by limiting K AOM to positive values. A different function K AOM (M, t^j is determined for each Monte 

Carlo trial. Because of interactions between the area oscillations model and the elastic modes of the parachute 
system (in the simulation due to the flexibility of the triple bridle), unrealistic values of the Mach number can be 
obtained at the parachute’s center of mass due to rapid motion of the parachute, resulting in nonphysical feedback 
with the area oscillations and Mach Efficiency Factor (section III.G) models. Thus, the Mach number used in the 
area oscillations model is that of the center of mass of the aeroshell, not that for the parachute center of mass. 

The equations used in the calculation of K AOM are given in the appendix. To illustrate the possible functional 
values that K AOM can take, two realizations are shown in figure 8 as a function time for a prescribed Mach number 
time history (also shown). The Mach number time history is that for the Viking BLDT AV-4 test flight, 11 starting 
near parachute inflation first peak force. Note that the calculated values of K AOM exhibits a random oscillatory 
behavior about a value of 1.0, its magnitude decreases with M , and becomes equal to one at M = 1.4 . 

H. Drag Coefficient and Mach Efficiency Factor 

A parachute drag coefficient, C Dq , model was created. This model considered the specific MSL geometry in its 
definition, including the effect of the longer parachute trailing distance (i.e., 10.32 •D AS , see figure 6). The 
reference areas for C Dq are the parachute nominal areas, either S 0 or (^o) Des ’ as appropriate (see sections III.B and 

III.F, and table 1). The parachute drag coefficient model is defined as functions of Mach number only. This model 
was created from an amalgamation of data including those from Viking wind tunnel tests, 17 Viking BLDT flight 
tests, 12,18 Phoenix Mars flight data, 16 and MSL-specific wind tunnel tests. Low , Nominal , and High values were 
defined and denoted as (C n ) , (C D ) , and (C n ) , respectively. The Low and High values define bounds 

V 0 / Low \ 0 'Nominal \ 0 / High 

of uncertainty surrounding (C n ) . These bounds are ±12.5 percent for M < 1, and increase asymptotically 

with Mach number for M > 1 from ±12.5 percent at M = 1 to ±25 percent at M = 2.6 . A graph of this C Dq model 
is shown in figure 9. Several features of the model can be observed. First, the model assumes that C Dq is 

independent of Mach number for M < 0.8 . Second, a transonic drag dip occurs near M = 0.95 . This behavior has 
been observed in wind tunnel test data, 17 and it is due to interaction between the aeroshell and parachute flow fields. 
Finally, after recovering from the transonic drag dip, the drag coefficient increases and peaks at a Mach number of 
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approximately 1.5, but then decreases again at higher Mach numbers. The C Dq model is only used directly in the 
simulation for the calculation of F PI Des (section III.F). However, it is used to define the Mach Efficiency Factor, 

MEF , model (see below) and is used to inform the static aerodynamic coefficients model (section III. I). 

The Mach Efficiency Factor, MEF , is calculated from the Nominal value of the drag coefficient model by the 
equation 


MEF = 


( c «.) 


Nominal 



Nominal @ M= 0.2 


(14) 


The Mach Efficiency Factor, as a function of the Mach number, is shown in figure 10. The purpose of the MEF is 
to serve as a scaling factor for the static aerodynamic coefficients to account for variations in their values as 
functions of M (section III.I). The behavior of MEF with respect to M is similar to that described above for C Dq . 

As mentioned in section III.G, there is an interaction between the area oscillations model, the MEF model, and 
the triple bridle elastic properties, that can yield unrealistic oscillatory values of the Mach number at the parachute’s 
center of mass due to rapid oscillations of the parachute. Thus, the values of the Mach number used with the drag 
coefficient and Mach Efficiency Factor models should be those of the center of mass of the aeroshell at all times. 


I. Aerodynamic Coefficients 

There are three static aerodynamic coefficients in the model created for the simulation: tangential (C r ) , normal 
(C N ) , and pitch (C mo ) . They are specified as functions of the total angle of attack, a PTot . The orientation of the 

forces and moments they represent are defined with respect to the total angle of attack plane as seen in figure 1 1 , 
where the force vectors specified by C T and C N are in the total angle of attack plane, and the moment vector 
specified by C mo is perpendicular to this plane. Because the parachute is axisymmetric, C T , C N , and C mQ are 
sufficient to describe the static aerodynamic coefficients of the parachute. As a set they can also be referred to as 
the total or axisymmetric static aerodynamic coefficients. The reference area used in their definition is S 0 . The 
reference length used for C mo is D 0 , and its reference point is the suspension lines confluence point. As defined 

here the suspension lines confluence point is the imaginary confluence point the physical suspension lines would 
reach if they were extended upstream until they met (see figure 6). In actual construction the suspension lines do not 
physically meet at this point since they are grouped into risers before they reach this point. If the normal force and 
the pitching moment are to be applied at a location other than this reference point (e.g., the parachute’s center of 
mass), the value of C mo needs to be recalculated for that point. 

The static aerodynamic coefficient model becomes active at suspension lines stretch and remains active 
thereafter (i.e., it is active for t > t LS ). As shown in figure 12, the static aerodynamic coefficients model is specified 
as a function of the parachute’s total angle of attack, at a Mach number of 0.2. For use at other Mach numbers the 
static aerodynamic coefficients need to be multiplied by the Mach Efficiency Factor defined in section III.H. Note 
that the Mach number to be used is that of the aeroshell (see discussion on this point in sections III.G and III.H). In 
the simulation the parachute total angle of attack and its dynamic pressure are calculated based on the airspeed 
components at the center of mass of the parachute (see table 3 and figure 6). During certain time intervals, and over 
certain values of Mach number, the static aerodynamic coefficients need to be multiplied by other constants: Kj 

during inflation (section III.F), K AOM during the time that area oscillations are possible (section III.G), by K MLE1 to 
account with interaction with the main landing engines (section III.K), and by K off to account for offloading 
(section III.L). 

Three curves are given for each of the static aerodynamic coefficients in figure 12: Nominal , Boundary 1, and 
Boundary 2. The Boundary curves define an estimated uncertainty surrounding the Nominal values. A uniform 
uncertainty distribution between the Nominal and Boundary 1 values, and between the Nominal and Boundary 2 
values, is used, with equal probability of being on either side of the Nominal value. The uncertainty boundaries are 
correlated. Note, for example, that the curve for Boundary 2 has the lowest value of the trim parachute total angle of 
attack (i.e., value of cx PTot for which C mo ), and presumably provides the best parachute stability. In agreement with 

the general trend that higher parachute stability is associated with a lower drag coefficient, it is the C T curve for 
Boundary 2 that has the lower values. Because of this correlation between curves, an interpolated set of static 
aerodynamic coefficient curves for a particular Monte Carlo trial must always be on the same side of the Nominal 
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curve and at the same relative distance. Such an interpolation can be accomplished by calculating an interpolation 
factor, K Int , defined by 

Kim = 2C/[0,l)-l (15) 


Thus, the range of K Int is -1 < K Int < 1 . Note that K Int is calculate only once for each Monte Carlo trial. The same 
value of K Int is used to interpolate all three static aerodynamic coefficients (i.e., {C T ) M , (C N ) lnt , and (C Wo )^ ? ). An 
interpolated set of curves for the static aerodynamic coefficients are then calculated from 


If K lnt > 0 then 


^t) m ~( C t) 

( C T)i m = ( C A 

{ C n)im = ( C A 

K )im ~ K ),\ 


Nominal ^ nt 


Nominal ^ nt 


+ K„ 


[(C T ) ~(C T ) 1 

\ 1 / Boundary 1 \ 1 / Nominal J 

[( C n) Boundary 1 " ) Nominal j 

\(c ) -(c) 1 

[V m oJ Boundary 1 \ m o J Nominal \ 


else 

,-K Inl [{C T ) 

r^[( c d. 

-K,„ 


— (C T ) 1 

Boundary 2 \ 1 ) Nominal J 

, y A c «) Nominal j 

lc ) -(c)! 

V m oJ Boundary 2 \ m o J Nominal J 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


endif 


The data used to create the static coefficient model was obtained during the test campaign described in 
reference 19. Judicious interpolation and extrapolation of these data was required to create this model. 

The static aerodynamic coefficient model is related to the drag coefficient model. At each of the three trim 
parachute total angles of attack where the pitching moment coefficient is zero, § and the normal force coefficient is 
also zero. Thus the associated values of the drag coefficients at the trim parachute total angles of attack are 
C Do = C T cosa P T ot • These three values of C Do calculated from the static aerodynamic coefficients for the 

Boundary 2, Nominal , and Boundary 1 curves were intended to be the same as the Low , Nominal , and High values 
of C Do , respectively, for M = 0.2 as determined from the drag coefficient model. An unintended small discrepancy 

of approximately 0.5 percent in absolute value exists between the Boundary l /High and Boundary 2/ Low values of 
C Dq at the trim parachute total angle of attack as determined per the above discussion. This was due to an oversight 

in updating models during development. However this discrepancy does not alter the fact that the drag coefficient 
and static coefficient models are related. 

No data was available to create models for the dynamic aerodynamic coefficients (e.g., C m ^ , C m . ). Thus, 

dynamic aerodynamic coefficients are not included as part of the parachute models and simulation. Doing so is the 
equivalent of assuming they are all zero. 


J. Apparent Mass 

For purposes of this paper, “apparent mass” is the term used to describe the mass of the fluid within and around 
the parachute that generates aerodynamic forces proportional to the parachute’s accelerations (both linear and 
rotational). The apparent mass model described here consists of a single point mass coincident with the physical 
center of mass of the parachute. This apparent mass, m App , is dependent on the freestream atmospheric density at 

the parachute’s center of mass, , and the parachute’s nominal diameter, D 0 . The apparent mass is calculated 
from 

m A PP = K AppJZp^D^ (22) 

where K App is the apparent mass coefficient, given by 

K App =T[ 0, 0.035, 0.070] (23) 


§ The trim parachute total angles of attack are 16.56° for the Boundary 1 curve, 8.90° for the Nominal curve, and 6.81° for the 
Boundary 2 curve. 
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The parameter T[0, 0.035, 0.070] is a random number with triangular distribution whose lower limit is 0, its mode 

is 0.035, and its upper limit is 0.070. A new value of T is determined for each Monte Carlo trial. Note that the use 
of T makes m App a perturbed parameter. The range of values for K App was selected after considering those used in 

reference 20, and the apparent mass discussion in course notes by S. Lingard. 11 In the simulation the apparent mass 
is subject to the force of gravity - an undesired effect since the apparent mass is only intended to generate inertial 
forces . To counteract the force of gravity on the apparent mass an opposite buoyancy force vector, F s , of 
magnitude 

Fg = m App§Mars ( 24 ) 

(where g Mars is the local acceleration of gravity) acts on the parachute’s center of mass. The apparent mass model 
becomes active in the simulation when the parachute reaches its first full inflation, and remains active thereafter 
(i.e., t>t FI ). 

The approach described here to model the apparent mass is highly simplified. See section V.C for 
recommendations regarding improvements to this model. 


K. Main Landing Engines Interaction 

The plumes of the main landing engines (MLE) can disturb the flow field experienced by the parachute. The 
Main Landing Engines Interaction (MLEI) model accounts for this interaction by defining a time-varying parameter, 
K MLE i > used to as a multiplier for the static aerodynamic coefficients. The value of K MLEI depends on: the 
location of the parachute with respect to the descent stage center of mass as shown in figure 13, the sum total thrust 
level of the MLE ( T MLE ), the dynamic pressure at the parachute, the nominal area of the parachute, the Nominal drag 
coefficient of the parachute, and a thrust scaling factor ( rj MLEI ). The MLEI model becomes active at the time MLE 
purging starts and remains active through the remainder of the EDL sequence. 

The volume within which the MLE interaction occurs is shown in figure 13. This volume is a cylinder of radius 
v mlei an d height h MLEI ( r MLEI = 80 m , h MLEI = 800 m ). One of the ends of this cylinder is centered at the parachute 
apex. The cylinder axis of symmetry is parallel to, and in the same direction as, the airspeed vector of the 
parachute’s center of mass, \ Poo . The MLE interaction will depend on whether the center of mass of the descent 

stage is inside or outside this cylinder. Two quantities describe the position of the descent stage center of mass with 
respect to the cylinder as shown in figure 13: r p DS and h P DS . These quantities can be calculated from: 


V P,DS 



(25) 


h 


P, DS 


V *R 


(26) 


where R p DS is the vector from the apex of the parachute to the center of mass of the descent stage as shown in 
figure 13. Note that r p DS is always positive, but h p DS can be either positive or negative. The value of K MLE1 is 
then calculated as follows: 

R f 3 , DS > r MLEI hp, DS > h ML El h p DS - < 0 thdl 


K MLEI ~ 1 

(descent stage center of mass is outside the cylinder) 
else 


(27) 




~ 1 V MLEI 


<?cA(Q>„) a 


(descent stage center of mass is inside the cylinder) 


endif 


(28) 


^ Lingard, S., “Aerodynamics 2 (Unsteady),” Lecture notes presented at the Parachute Systems Technology Short Course , U.S. 
Army Yuma Proving Ground, May 12-16, 2008. 
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If K mlei < 0.02 then K MLEI = 0.02 endif (29) 

(limits the minimum value of K MLE1 to 0.02) 

where rj MLEI is a scaling factor controlling the magnitude of the MLE interaction. The value of rj MLEI was set at 
0.3075. Note that a minimum value for K MLEI has been enforced. The possible range for K MLE1 is 0.02 < K MLE1 < 1 . 
The static aerodynamic coefficients of the parachute are multiplied by K MLE1 before using them to calculate the 
aerodynamic forces and moments on the parachute. 

The MLEI model was included in the simulation to address concerns regarding backshell re-contact with the 
separated descent stage and/or rover, either shortly after separation during MLE throttle-up, or long after separation. 
The model is derived from conservation of momentum principles. However, not enough data are available to assess 
the validity of its physics and to properly tune the model (i.e., for selecting values of rj MLEI , h MLEI , and r MLEI ). In 
the absence of sufficient knowledge regarding this model, the simulation approach was to bound the problem by 
incorporating a large interaction cylinder, and exercising the simulation both with and without the MLEI model to 
assess the short- and long-term re-contact risks. 

L. Offloading 

Experimental data from reference 2 1 indicates that a sudden offloading of the parachute due to a step reduction 
in payload mass (i.e., lander and/or descent stage separation) results in a temporary reduction in the drag area of the 
parachute and backshell combination. Analyses indicated that this reduction in drag area could not be accounted for 
by the modeled parachute apparent mass (section III.J). Thus, to account for this temporary reduction in drag area, 
the offloading model described here was created and implemented in the simulation. This model was created by 
reanalyzing the data in reference 21 for MSL. Thus this model is not general; application to other spacecraft would 
require reanalysis of the source data in reference 2 1 . 

The parachute offloading model takes the form of a time-dependent nondimensional parameter, K off , used as a 

multiplier for the parachute and backshell static aerodynamic coefficients. This multiplier becomes active at the 
time of descent stage separation from the backshell, t = t Sep DS . The offloading multiplier, K off , is calculated from 

^Ojf = 1 “ (l “ ) eX P [“ (f ~ hep, DS ) 1*0]! ] 0 0) 

where £ Min =0.793 and r off =1.965 s . Both t and t Sep DS must be entered in seconds in equation (30). A plot of 
K off vs. t-t Sep DS for 0 s ^t-t Sep DS < 6 s is shown in figure 14. Note that within three seconds of t Sep DS , K off is 
greater than 0.95, and within six seconds it is greater than 0.99. 

IV. Models Creation, Implementation, Checkout, and Configuration Control 

During the process of developing the models described herein, several lessons were learned regarding their 
creation, implementation, checkout, and configuration control. Some of these lessons are listed below. 

• The responses (i.e., results) required from the simulation need to be clearly understood. Models created for 
the simulation should be focused on obtaining accurate and/or conservative values of these responses. 
Models that do not support these responses in a meaningful way should be considered for deletion or 
simplification. 

• Some models are conservative with respect to some simulation responses, but unconservative with respect to 
others. For example intentionally overestimating the drag coefficient may be conservative with regards to 
opening forces, but unconservative with regards to terminal descent. In such situations a risk assessment may 
need to be conducted to balance the risks throughout the system. Another possibility is to try and separate 
aspects of the model (or perhaps create separate models) so that conservatism can be tuned into the model(s) 
to address the risk posture adopted by the project regarding various responses. To continue with the previous 
example, the drag coefficient may be overestimated in the Mach number range where inflation is expected to 
occur, but underestimated at the lower Mach numbers associated with terminal descent. 

• Models should be clearly documented. This documentation should be sufficient for the model to be 
implemented with minimal interaction with the person(s) that created and documented the models. The 
model documentation should have a revision history so that changes made to the models can be traced. 
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• After the models are implemented into the simulation, the person(s) that created the models and the person(s) 
that implemented them in the simulation should work together to verify that the modes have been 
implemented and function as desired. 

• A system for configuration control of the models and their data should be established and used. For models 
that appear in multiple documents, the governing (ruling) source for the models and/or data should be clearly 
established. 


V. Model Improvements 

Based on the experience of creating the models described above, implementing them in the simulation, and/or 
comparing the model predictions to the reconstructed values from the MSL flight data, 4 several improvements to the 
parachute system models can be suggested. These improvements are described below. 

A. Deployment Time 

The present deployment time model defined by equations (3) to (5) in section III.F is empirical and based on 
data from previous missions, in particular the Viking BLDT flight tests 9 " 11 and Mars landers, 13 and the Phoenix Mars 
lander. 14 This was appropriate because of the similarities in the parachute system and operating conditions among 
these spacecraft. Reconstruction of the MSL flight data 4 indicates that the model over-predicted the deployment 
time. It is possible to create a model that explicitly took into account the mortar system performance, aeroshell 
deceleration at mortar fire, and other relevant variables to create a physics-based model that could predict the 
deployment time more accurately. 

B. Wind Angles Calculation 

The parachute total angle of attack used as the independent variable in the static aerodynamic coefficient model 
was calculated in the simulation at the parachute’s center of mass. It is reasonable to assume that the parachute’s 
aerodynamic forces are closely related to the wind angles somewhere near the canopy. However, the parachute’s 
center of mass - where the wind angles are currently being calculated - is 7.8 m upstream of the parachute’s skirt 
(see figure 6). In addition, the parachute’s center of mass can change position due to changes in materials, without a 
concurrent change in the inflated shape of the parachute. Thus the parachute wind angles should be calculated at a 
fixed geometric location near the canopy, perhaps at the skirt. 

The difference between the wind angles calculated at the parachute’s center of mass and the skirt was estimated 
to be less than one degree for MSL. Thus, in the simulation it was decided not to change the location where the 
wind angles were calculated. However this is not a general conclusion, and the suggested improvement would be 
relatively simple to implement (it was not done so for the simulation to avoid having to revalidate the already tuned 
simulation late in its development, and to prevent unintended consequences). 

C. Apparent Mass 

The apparent mass model that was used in the simulation is highly simplified. The apparent mass cannot be 
accounted for with a single point mass - it has different values depending on the axis, and has both mass and mass 
moment of inertia components. It is recommended that a more sophisticated apparent mass model be incorporated 
into the simulation for future use, perhaps using the approach suggested by S. Lingard (see footnote U in 
section III.J). 


VI. Concluding Remarks 

Parachute-related models (see figure 4) were incorporated into the MSL POST2 EDL simulation. These models 
were kept as simple as possible, considering the overall objectives of the simulation (section I). This paper has 
described the models to the extent necessary to understand how they work and some of their limitations. In a 
companion paper, 4 results from a parachute performance reconstruction using Mars flight data for MSL are 
presented, and compared to some of the parachute models and simulation discussed here. This comparison indicates 
that the parachute models and simulation provided a reasonable prediction of the parachute system flight behavior 
on Mars. 
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Appendix: Area Oscillations Model Equations 

The equations used in the calculations to obtain the area oscillations multiplier, K AOM , as described in section 
III.G are given below. The data used to tune this model were obtained from the Viking BLDT AV-4 flight test 11 and 
the Phoenix Mars flight. 14 " 16 The variable M is the Mach number. The variable M FI is the Mach number at full 

inflation. Because of interactions between the area oscillations model and the elastic modes of the parachute system, 
unrealistic values of the Mach number can be obtained at the parachute’s center of mass due to rapid motion of the 
parachute. Thus, the values of M FI and M used in the area oscillations model should be those for the center of 
mass of the aeroshell, not those calculated at the parachute’s center of mass. The model should not be used for 
Mach numbers higher than M Init . 
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Constants used in the area oscillations model. 


Constant 

Value 

n AOM 

3 

k HFi 

0.20 

k HF2 

0.15 

Min 

2;r/12~ 0.524 rad/s 


2.057 

M final 

1.40 

E Max 

0.59 

T High 

1.08 

^Low 

1.00 


A = 


In 2 


M 


Init 


-M 


Final 


K High =15.9474-6.5617M FI 

K Low = 2.43032 + 0.787402M f/ 

Rj = independent C/[0,l) for every i = 1 to 6 
values for R { ’s calculated once per Monte Carlo trial 


*s = 



K t = R s R j for i = 1 to 6 

co i = ;r{f/[0,l) + (/ - 1)} + co Mn for every i = 1 to 6 
is independent for every i 

values for co i ’s calculated once per Monte Carlo trial 
units of co i : rad/s 

= 2tt£/[ 0,1) for every i = 1 to 6 
is independent for every i 

values for ip t ’s calculated once per Monte Carlo trial 
units of ip i : rad 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


(A8) 
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co hf 1 -^{7t/[0,l) + 6} + w Mn 
^HF 2 = ft [o, l) + 1 3 j + CO Min 

£/[0,l) is independent for co HF and (o HF 
values for co HF and cd hf ^ calculated once per Monte Carlo trial 
units of co np and co HF : rad/s 


t shift =t_t Fi where t is the lowest value of t > t FI for which: 


J^.sin^[^(r-f F/ ) + ^] =0 


and its derivative 



n AOM 2 w i K i sin"' 10 *'" 1 [«, {t-t F ,) + n>] cos [m, (t - t FI ) + xp t ] < 0 


(A9) # 


(A10) 


(All) 


At this point all constants have been determined, and the calculations proceed as a function of time for t > t FI . 


6 


S « (0 “ 2 ^ Sin ” A0M M ? + 1 shift - { FI ) ■ + n>i ] 

i= 1 

(A 12) 

^6 (0 ^ hf ] S ^ n AOM \ C 0 hf 1 tpi )] K HFi sin A0M [co HFi (t A/)] (A 13) 

If T x (£) > 0 then 


T 2 (0 = 1 + K mgh f (0{exp[A(M - - 1 

[ (A 14) 

f (0 = 1 + E Max (l - expj-r^ [T 2 ( t ) - 1]}) 

(A 15) 

else 


T 2 (0 = 1 ■ + KusJi (0{ ex P [ A (M - M Final )] - 1 J 

(A 16) 

7’ 3 (0 = exp{T to [7’ 2 (0-l]} 

(A 17) 


endif 

perform the above conditional evaluation and appropriate calculations at every time step 

If M > M Final then K AOM ( t ) = fit) else K A0M (t) = 1 endif (A18) 

perform the above conditional evaluation and appropriate assignment at every time step 


# All quantities with the subscript “HF” are related to the high-frequency terms shown in equation (A 12). 
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Table 1 . Key Geometric Parameters of the MSL Parachute System 


Item 

Symbol 

Value 

Comments 

Nominal (reference) parachute 
diameter 

D 0 

21.348 m + (0.2 m-U[0,l) - 0.1 m) 

The unperturbed value 
(21.348 m) is that of the flight 
unit parachute. Random 
portion is for use in Monte 
Carlo trials. 

Nominal (reference) parachute 
area 

So 

S 0 = jtD 0 / A 

Depends on perturbed value of 
D 0 . Unperturbed value: 

S 0 = 357.94 m 2 . 

Nominal (reference) parachute 
diameter used in the 
calculation of F fi Des only 

(D<> ) Des 

21.50 m 

No perturbation. 

Nominal (reference) parachute 
area used in the calculation of 

F FI,Des ° nl y 


363.05 m 

Calculated using the equation: 
(So) Des = x(D 0 ) 2 j4. 

Unstretched (as-built) length 
of the triple bridle legs (each) 

LtBL 

1.385 m 

See figure 6. 

Distance between the triple 
bridle confluence point and 
attachment points on the 
aero shell (unstretched, as-built 
triple bridle legs) 

H tb 

1.349 m 

See figure 6. 

Single riser length 

Lr 

7.546 m 

See figure 6. 


Notes: 

• The quantity U[0,1) is a random variable with uniform distribution in the interval between 0 (inclusive) and 1 
(not inclusive). Unless specifically called for, the value of U[0,1) is not reused. 
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Table 2. As-built key geometric parameters of the MSL parachute. 


Parameter 

Value 

Nominal (reference) diameter, D 0 (m) 

21.348 

Disk diameter, D d (m) 

15.447 

[D d /D o] 

[0.724] 

Band diameter, D B (m) 

15.608 

[D b /D 0 ] 

[0.731] 

Vent diameter (unconstrained), D Vu (m) 

1.576 

[DJD 0 ] 

[0.0738] 

Vent diameter (constrained), D Vc (m) 

1.524 

[D Vc /D 0 ] 

[0.0714] 

Gap height, H G (m) 

0.904 

[H g /D o] 

[0.0424] 

Band height, H B (m) 

2.580 

[H b /D 0 ] 

[0.1209] 

Suspension lines length, L s (m) 

36.578 

[L s /D 0 ] 

[1.713] 

Projected diameter, D p (m) 

14.944 

[D P /D 0 ] 

[0.700] 

Nominal (reference) area, S 0 (m 2 ) 

357.94 

Disk area, S D (m 2 ) 

187.40 

[S D /S 0 ] 

[0.524] 

Gap area, S G (m 2 ) 

44.01 

[So/So] 

[0.123] 

Band area, S B (m 2 ) 

126.52 

[S B /S 0 ] 

[0.353] 

Vent area (unconstrained), S Vu (m 2 ) 

1.95 

[SJS (] ] 

[0.0055] 

Vent area (constrained), S Vc (m 2 ) 

1.82 

[AAJ 

[0.0051] 

Geometric porosity (unconstrained vent), X gu (%) 

12.84 

Geometric porosity (constrained vent), A gc (%) 

12.80 

Number of gores (and suspension lines), N gores 

80 



Parachute constructed shape. 
Figure not to scale. 


Notes: 

• Dimensions in this table are, in general, not known to the level accuracy implied by the numbers above. They are 
presented here with extra significant figures to make them internally consistent. 

• Most of the data in this table was provided by A1 Witkowski of Pioneer Aerospace. 

• Constraint on vent diameter is due to the vent cords. 

• As defined here the suspension lines length is the distance from the parachute skirt to the imaginary confluence point the 
physical suspension lines would reach if they were extended. In actual construction the suspension lines length as defined 
here consists of both physical suspension lines and risers that collect sets of suspension lines into fewer load paths. 

• The projected diameter D p and the ratio D p /D 0 are estimates. 
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Table 3. Nominal Inflated Parachute Mass Properties. 


Property 

Symbol 

Value 

Mass 

m p 

46.826 kg 

Center of mass along the axis of symmetry, downstream from the 
suspension lines confluence point 

( X chute) CoM 

-27.516 m 

Mass moment of inertia about the parachute’s axis of symmetry 

( I xx ) CoM 

1,140.1 kg»m 2 

Mass moments of inertia about axes perpendicular to the 
parachute’s axis of symmetry and to each other, with origin at the 
parachute’s center of mass 

yy)coM an< ^ ( I zz ) CoM 

10,477.0 kg*m 2 


Notes: 

• The mass of the parachute stated here, m p , includes the deployment bag, canopy, suspension lines, suspension 
lines risers, and half the single riser. 

• The parachute axis of symmetry is x chute , with origin at the suspension lines confluence point, positive pointing 

upstream. As defined here the suspension lines confluence point is the imaginary confluence point the 
physical suspension lines would reach if they were extended upstream until they met. In actual construction 
the suspension lines do not physically meet at this point since they are grouped into risers before they reach 
this point. 

• The axes y chute and z chute are perpendicular to the parachute axis of symmetry, x chute , and to each other. 

• The parachute geometric model used to derive the mass properties in this table is axisymmetric. Thus, the 
parachute’s center of mass is on its axis of symmetry, the mass products of inertia perpendicular to the axis of 
symmetry are identical, and the mass products of inertia are zero. 
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Figure 1. MSL parachute in operation on Mars photographed by the HiRISE camera on the Mars Reconnaissance 
Orbiter. Image credit: NASA / JPL-Caltech / Univ. of Arizona. 



Figure 2. MSL parachute being tested at the National Full-Scale Aerodynamics Complex wind tunnel located at the 
NASA Ames Research Center in Moffett Field, California. The anti-inversion netting seen attached to the skirt of 
the canopy was only present for subsonic testing - it was not present on the parachute flown on Mars. Image credit: 
NASA / JPL-Caltech. 
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Figure 3. MSL EDL concept of operation from just prior to parachute deployment. Graphics source: JPL-Caltech. 
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Parachute Models 


EDL Events 



Start SUFR maneuver 
RCS damping enabled 

RCS damping inhibited 
Mortar fire 

Start parachute deployment, , 
Parachute suspension lines st 

Parachute 1 st peak force 
M= 1.4 

RCS damping enabled 

RCS damping inhibited 
Heatshield release, 0.58 

RCS damping enabled 

MLE purging initiation 

RCS damping inhibited 
Descent stage separation 


Figure 4. Relationship between the parachute models and EDL events. A closed circle at the beginning or end of a 
model vertical line indicates that the model is active at that event. A closed circle at the end of a horizontal dashed 
line indicates a specific time. A single arrow on the end of a model vertical line indicates that the model is active up 
to, but not including that event. A double arrow on the end of a model vertical line indicates that the model is active 
until landing. Mass properties (A) are those retained by the aeroshell; Mass properties (B) are those associated with the 
parachute. See the specific model sections for more details. 
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Figure 6. Parachute system dimensions. All linear dimensions in meters. Note that the triple bridle legs attach to 
the aeroshell inside its outer mold line. Aeroshell graphic by Karl T. Edquist, NASA Langley Research Center. 
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K:=C X 


t 


Figure 7. Graphical representation of the inflation multiplier Kj . 



Figure 8. Two examples of the area oscillation multiplier, K AOM , vs. time since mortar fire, t — t MF , for an assumed 
Mach number, M , time history (Viking BLDT AV-4 test flight 11 ). 
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Figure 9. Drag coefficient, C Dq , model. 
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Figure 11. Parachute static aerodynamic coefficients. 





Figure 12. Parachute static aerodynamic coefficients model. The trim parachute total angles of attack are 16.56° 
for the Boundary 1 curve, 8.90° for the Nominal curve, and 6.81° for the Boundary 2 curve. 


26 

American Institute of Aeronautics and Astronautics 



Figure 13. Geometry for the Main Landing Engines Interaction (MLEI) model. As shown here, h p DS > 0 . Figure 
not to scale. 
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